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Abstract 

In this paper, we study optimal transportation problems for multifractal random 
measures. Since these measures are much less regular than optimal transportation 
theory requires, we introduce a new notion of transportation which is intuitively 
some kind of multistep transportation. Applications are given for construction of 
multifractal random changes of times and to the existence of random metrics, the 
volume forms of which coincide with the multifractal random measures. This study 
is motivated by recent problems in the KPZ context. 
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1 Introduction 

On the Borelian subsets of M m , consider a measure M formally defined by 

(1.1) M(A)= [ etXto-^xWdx, 

J A 

where 7 > is a parameter, (X(x)) x£ ^m is a Gaussian distribution with covariance 
function given by 

(1.2) Cov(X(x),X(y)) =ln + -^—+g(\y-x\), 

\y - x\ 

and g is a continuous bounded function. Actually, whatever the function g is, the part 
that really matters in (1.2) is the logarithm part. Such measures are called log-normal 
multifractal random measures and were first rigorously defined in [4]. 
The above situation can be generalized to the situation 

(1.3) M(A)= [ e^ x) dx 

J A 
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where the process u is a suitable Levy distribution: the resulting measures are called 
log-infinitely divisible multifractal random measures, MRM for short (see subsection 2.1 
for a reminder of the construction). Such measures exhibits interesting properties like 
stationarity, isotropy, long-range dependence, fat tail distribtions and, because of their 
log part in (1.2) (or a suitable generalization for Levy distributions), possess a remarkable 
scaling property, the so-called stochastic scale invariance: 

(1.4) (M(XA)) AcmT) lawa ^° (\^M(A)) AcmT) , 

where Q\ is an infinitely divisible random variable, independent of (M(A)) AcB ( 0T y 
When M satisfies the above relation (1.4) with = instead of ~, we will say that M 
satisfies the exact stochastic scale invariance property (or M is ESSI for short). 

The purpose of this paper is to investigate (optimal) transportation problems asso- 
ciated to these measures. A transport map between two probability measures /x, v is a 
map that pushes fi forward to v. The transport map is said to be optimal if it real- 
izes the infimum of a cost functional among all the possible transport maps. For usual 
cost functionals, existence and uniqueness of an optimal transport map are strongly con- 
nected to the regularity of the measures fi, v. Concerning MRM, their regularity is much 
weaker than that required in optimal transportation theory. So we give new notions of 
transportation that can be applied to MRM. This study is highly motivated by the con- 
struction of multifractal random changes of time and the construction of metric spaces 
the volume form of which is given by the MRM, as explained in subsection 3. The latter 
construction is related to a recent problem in the KPZ context raised by K. Izyurov at 
the conference "Les Diablerets". 

2 Background and main results 

In this section we first give the basic background in order to state rigorously the main 
results of the paper. Since the function g in (1.2) does not play a part in what follows, 
we focus on the case where the measure M satifies the exact scale invariance property. 

2.1 Reminder of the construction of ESSI MRM 

We present below the generalization of (1.1) to the situation where X is a Levy distri- 
bution. For further details, the reader is referred to [8]. To characterize such a Levy 
distribution, we consider the characteristic function of an infinitely divisible random 
variable Z, which can be written as K[e iqZ ] = where (Levy-Khintchine's formula) 

(p(q) = ibq - \a 2 q 2 + f (e iqx - 1 - iqsin(x)) v(dx) 

and u(dx) is a so-called Levy measure satisfying J R » min(l, x 2 ) v{dx) < +oo. We also 
introduce the Laplace exponent ip of Z by ip(q) = ip(—iq) for each q such that both terms 
of the equality make sense, and we assume that = (renormalization condition), 
-0(2) < +oo and Jr_ 1 ^ \x\ v{dx) < +oo (sufficient conditions for existence of MRM). 

Now we define the process u of (1.3). We remind that this process has to be stationary, 
isotrop and suitably scale invariant. To ensure isotropy, we introduce the unitary group 
G of M m , that is 

G = {M € M m (R); MM 1 = I}, 



2 



and H its unique right translation invariant Haar measure with mass 1 defined on the 
Borel o"-algebra B(G). For stationarity and scale invariance, we consider the half-space 

S = {(t,y);t 6 1,1/eR;}, 

with which we associate the measure (on the Borel er-algebra B(S)) 

6(dt,dy) = y~ 2 dtdy. 

Then we consider an independently scattered infinitely divisible random measure \i as- 
sociated to (<p,H ® 6) and distributed on G x 5. In particular, we have E(e l9/1 ^) = 

We combine the previous ingredients in a geometrical way to define uj. Given T, let 
us define the function / : M + — > M. by 

(l)= (l, iU<T 
n ' \ T if I > T 

The cone-like subset A[(t) of S is defined by 

Mt) = {(s, y)eS;y>l, -f(y)/2 <s-t< f(y)/2}. 

For any x G M m and g G G, we denote by xf the first coordinate of the vector gx. The 
cone product Ci(x) is then defined as 

d{x) = {(g,t,y) GGx S;(t,y) € Mxf)}, 

and the process (0 < I < T) by ui(x) = fi(Ci(x)) for x G W 71 . 

The Radon measure M is then defined as the almost sure limit (in the sense of weak 
convergence of Radon measures) by 

M(A) = lim MA A) = lim / e^ x) dx 

for any Lebesgue measurable subset A C R m . The convergence is ensured by the fact that 
the family (M;(A));>o is a right-continuous positive martingale. The structure exponent 
of M is defined by 

> 0, ((q) =dq- 

for all q such that the right-hand side makes sense. The measure M is different from 
if and only if there exists e > such that ("(1 + e) > m, (or equivalently ip'(l) < m). In 
that case, we have: 

Theorem 2.1. The measure M is stationary, isotropic and satisfies the exact stochas- 
tic scale invariance property: for any A g]0, 1], 

(M(XA)) AcB{0tT) la = (X m e^M(A)) AcB(0tT) , 

where Q\ is an infinitely divisible random variable, independent of (M(^4)) 7 4 CB ( ,T) > the 
law of which is characterized by: 

E[e iqQx } = 
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Furthermore, the support of such measures is full in the sense that 

Proposition 2.2. If the measure M is non-degenerate, that is < m, we have 

Supp(M) = M. m . Consequently, every Borelian subset ofW" 1 with null M-measure has 
its complement dense in R m for the Euclidian distance. 

Proof. For a given ball B(x,r) C M. m , the event {M(B(x,r)) > 0} is measurable with 
respect to the asymptotic sigma algebra generated by (w/)z>o- By the 0-1 law, it has 
probability or 1. Because of the uniform integrability of the martingale (Mj(j4))/>o 
for each Borelian subset A with finite Lebesgue measure (denoted by X(A)), we have 
E[M(B(x,r))} = X(B(x,r)) > 0. We deduce F(M(B(x,r)) > 0) = 1. Hence, P a.s., 
M(B(x,r)) > for all the balls B(x,r) with rational centers and radii. □ 

2.2 Optimal transport for MRM 

The reader may consult Appendix A for a brief reminder about optimal transport theory. 
In particular, we adopt the notations used in Appendix A. We denote by Br the closed 
ball of M m centered at with radius R and by Cr its Lebesgue measure. We define Xr as 
the renormalized Lebesgue measure on Br with mass 1 and the renormalized probability 
measure on Br 

mdx) = M(k) M(£te >' 

Originally, our problem was to construct a nice (as much as possible) mapping that 
pushes the MRM M forward to the Lebesgue measure. Of course, there are several 
possibilities to achieve that construction. Among all the possibilities, we focus on optimal 
transportation theory because it features many interesting qualities: measurability with 
respect to the randomness, isotropy, regularity,... Basic results ensure that there exists an 
optimal transport map pushing M forward to the Lebesgue measure Xr (for a quadratic 
cost function) provided that the measure M does not give mass to small sets. Though 
that condition is usually not satisfied by MRMs, this is true when the intermittency 
parameter ^'(1) is small: 

Theorem 2.3. When ip'(l) < 1, the measure M does not give mass to small sets. Hence 
there is a unique optimal transport map (for the Euclidian quadratic cost) that pushes 
the renormalized probability measure M forward to the renormalized Lebesgue measure. 

We can thus apply the classical transportation theory (see Theorem A. 3). There exist 
two optimal transport maps 

X ■ Supp(x) -> Supp(r) and T : Supp(r) -> Supp(x), 

(the supports of which are Borelian and contained in Br) that respectively push the 
Lebesgue measure Xr forward to M and vice-versa, meaning 

X# X R = M and T # M = Xr. 

They are both unique, gradients of convex functions and bijections, inverse from each 
other. They are respectively Xr and M-almost surely defined, meaning M(Supp(r)) = 1 
and A#(Supp(x)) = 1, so that both supports Supp(r), Supp(x) are dense in Br for the 
Euclidian distance. 
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(a) 7 2 ~ 0.1 



(b) 7 2 ~ 1.5 




(c) 7 2 ~ 2 (d) 7 2 ~ 3.9 

Figure 1: Simulations of the density of a log-normal MRM with various intermittency parameters 

2 2 

7 appearing in (1.1). The structure exponent matches £(g) = (2 + ^-)g — "2"<Z 2 - I n dimension 2, 
the measure is non-degenerate provided that < j 2 < 4. The last two simulations are colored 
with a logarithmic intensity scale. 
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We equip Br with the Riemannian metric 



M(Br) 2 

V* e B R ,Vu,v e R m ,g e x (u,v) = { R) (u,v) 

where (•, ■) denotes the usual inner product on R m . In that way, (Br, g e ) is a Riemannian 
space in which the volume form matches M(Br) x Ar(cIx) and the geodesic distance d e 
is given by the Euclidian distance on Br up to a multiplicative constant: 

d e{x,y) = M J^ R \ x - y\. 

On the support Supp(r) of T, we can define the distance dr by 

Vx,y 6 Supp(r), d r (x,y) = d e (T(x),T(y)). 

Hence there is an isometry, namely T, between the metric-measure space (Supp(r), dr, M) 
and the metric-measure space (Supp(x), d e , M(Br) x Xr). Since the closure of Supp(x) 
with respect to the Euclidian distance is equal to Br, the completion of the metric 
space (Supp(r), dr), denoted by (C, dr), is isometric to (BR,d e ). That isometry, which 
coincides with T on Supp(r) is still denoted by T and its inverse, which coincides with \ 
on Supp(x), is still denoted by \- Obviously, the metric space (C, dr) is compact. 

Furthermore, since Supp(r) is a Borelian subset of Br as well as a Borelian subset 
of C (for the respective topologies of Br and C) and since M(Supp(T)) = M(Br), the 
measure M can be extended to the whole of the Borelian subsets of C by prescribing: 

for any Borelian subset A of C, M(A) = M(A n Supp(r)). 

Via pullback, the space C inherits the structure of Riemannian manifold (smooth, 
complete, connected, m-dimensional manifold equipped with a smooth metric tensor). 
The (only) chart is given by V : C —¥ Br. We summarize below what we have proved as 
well as the properties inherited from the pullback metric: 

Theorem 2.4. Ifip'il) < 1, we can find a compact Riemannian manifold (C,g) and a 
Borelian subset B of Br such that: 

1. B is dense in Br for the Euclidian distance and has full M -measure, namely 
M(B R \ B) = 0, 

2. C is the completion of B with respect to the geodesic distance on C , 

3. the volume form on C coincides with the measure M on B, 

4- in the system of local coordinates given by the chart T, the Riemannian metric 
tensor on C reads 

M{Br) 2 

g = d(uj)(dx\ + • • • + dx^) with d(ui) = — — ^ • 
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The strongly intermittent case 1 < tp'(l) < m 

Now we consider a log-infinitely divisible random measure M satisfying 1 < if)'(l) < m. 
The measure M gives mass to small sets so that we cannot use classical theorems of 
optimal transport theory if we consider the ambiant space Br equipped with its Euclidian 
structure. It thus seems hopeless to solve the problem: 

(2.1) Find the mapping <p realizing the infimum: inf / \(p{x) — x\ 2 M(dx). 

f-B R -^B R J Br 

<P#m= X R 

Our strategy is the following: can we find a pair (M',T'), where M' is a probability 
measure and T' is an optimal (for the Euclidian quadratic cost) transport map pushing 
M' forward to Xr, such that the following Monge type optimization problem possesses a 
(unique) solution? 

(2.2) Find the mapping ip realizing the infimum: inf / \(p{x) — T' (x)\ 2 M (dx) . 

f-B R -+B R J Br 

It turns out that the above problem is a mathematical formulation of the following 
intuitive argument: if we cannot find an optimal transport pushing M forward to Xr, can 
we find an intermediate measure M' and optimal transports T, T' respectively pushing 
M forward to M' and M' forward to Xr? If the answer is positive, then the composition 
T'oT pushes M forward to Xr. Though the composition T' oT is in general not optimal 
for the Euclidian quadratic cost, it is the composition of two gradients of convex functions, 
which is not bad in terms of regularity. So we have formalized some kind of two-step 
optimal transport. And more generally, if we cannot find a two-step optimal transport, 
is it possible to find a n-step optimal transport between M and Xr, that is a composition 
of n gradients of convex functions? 

The reader may have the following objection: from classical theorems, existence of a 
(unique) optimal transport map pushing M forward to M' or Xr does not depend on the 
target measure (M 1 or A_r) but only on M through the fact that M does or does not give 
mass to small sets. So, a priori, two-step optimal transports may be as difficult to exhibit 
as optimal transports. Our idea lies in the fact that Problem (2.2) can be reformulated 
in a very simple way if we equip the ball Br with an appropriate Riemannian structure. 
Indeed, if we change the unknown in (2.2), we get the following equivalent problem 

(2.3) 

Find the mapping if) realizing the infimum: inf / \T'(if>(x)) — T' (x)\ 2 M (dx) . 

tP:B r ^B r J b 
^M= M ' 

It turns out that we can equip the ball Br with a Riemannian structure, the distance 
of which matches d(x,y) = \T'(x) — T'(y)\ and the volume form of which matches M' 
(up to a multiplicative constant). Problem (2.3) thus reduces to a classical problem of 
optimal transportation theory on smooth Riemaniannian manifolds: for such a mapping 
to exist, it is mainly sufficient that M does not give mass to the small sets associated 
to the distance d. That is the main constraint when choosing the measure M': it must 
be the volume form associated to a Riemannian metric, the small sets of which are not 
charged by M. Of course, the argument generalizes to n-step optimal transports. 
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The main difficulty thus lies in choosing the number n of steps and the intermedi- 
ate measures. The crucial point is the following. Given n > 1, M can be seen as the 
composition of n multiplicative chaos (see subsection 4.2): we can find n independent in- 
dependently scattered log-infinitely divisible random measures ■ ■ ■ , // n ^ associated to 
(<p/n, 9) (see subsection 2.1). The corresponding processes ui associated to . . . , 
are respectively denoted by oj\ , , . . . , uj\ . . We define recursively for k < n: 

M {0) (dx) = dx and M {k) (dx) = lim e^^M^'^idx), 

where the limits have to be understood in the sense of weak convergence of Radon 
measures. Then both measures M and have the same law so that we assume, in 

what follows, that M and coincide. That procedure allows to see each measure 

chaos with respect to M^ k ^ with a reduced intermittency parameter ip'(l)/n. 
If we can equip the ball Br with a Riemannian metric g( k ~^ the volume form of which 
coincides with M^ k ~ l \ it turns out that does not give mass to the g( fc_1 )-small 

sets provided that the intermittency parameter tfj'(l)/n is small enough. So it suffices to 
choose n small enough. In consequence there is a unique optimal transport map (w.r.t. to 
the quadratic cost function associated to the metric g^ k ~^) that pushes the renormalized 

measure M forward to the renormalized measure M . And so on for the different 
values of k < n. Thus we claim: 

Theorem 2.5. If 1 < tp' (1) < d, we can find n > 1 and n gradients of convex functions 
. . . , T( n ) such that \/k = 1, . . . , n, the mapping = o ■ • • o pushes the 

measure M forward to the Lebesgue measure and minimizes the quantity 

inf / \T(x) -ip {k - l) {x)\ 2 M {k \dx). 



T:B R ->B R j Br 
T #I jC=)- A « 

As a corollary, we get 

Theorem 2.6. We can find a compact Riemannian manifold (C,g) and a Borelian subset 
B of Br such that: 

1. B is dense in Br for the Euclidian distance and has full M -measure, namely 
M(B R \B) = 0, 

2. C is the completion of B with respect to the geodesic distance on C. 

3. the volume form on C coincides with the measure M on B 

4- in the system of local coordinates given by the chart ip( n \ the Riemannian metric 
tensor on C reads 

M (Br) 2 

g = d(uj)(dx\ + • • • + dx^) with d(uj) = — — ^ • 

Remark 2.7. In Theorems 2.4 and 2.5, if R < T, by scale invariance, the random 
variable d has the same law as e % where Q.r is an infinitely divisible random 

variables the law of which is characterized by^,[e Q t] = (^) So the radius R of the 

ball Br influences the random metric through a log-infinitely divisible random variable, 
which turns out to be log-normal for log-normal MRM. 
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3 Applications 



3.1 Application to the KPZ framework 

A special case of (1.1) in dimension 2 has recently received much attention. When X is a 
Gaussian Free Field (GFF), that is a Gaussian process with covariance function given by 
the Green function of the Laplacian, the measure M is called quantum measure. That 
situation is motivated by Liouville quantum gravity [3]. Much effort has recently been 
made to understand the geometry of these measures. The first important step was to 
prove the so-called KPZ formula [2, 3, 7]. Roughly speaking, the KPZ formula gives the 
correspondence between the Haudorff dimension of a set as seen by the Lebesgue measure 
and the the Haudorff dimension of this set as seen by the measure M. More precisely, 
for a given compact set E C Br and a mesure u, the s-dimensional Hausdorff measure 
of E wrt to v is the quantity: 

(3.1) 

H s (u, E) = lim HI(E, u), where 
5 — >-0 

H$(u, E) = inf{J^ y{B n ) s/2 ; E C |J B n , B n measurable set with < M{B n ) < 5}. 

n 

Then we define the Hausdorff dimension dim^(£') of the set E wrt to the measure v: 

(3.2) dim^(£) = inf{s > 0; H s {u, E) = 0} = sup{s > 0; H s (u, E) = +oo}. 

When v is given by the Lebesgue measure, the corresponding Hausdorff measures and 
dimensions will be called Euclidian. The KPZ formula asserts that the Euclidian Haus- 
dorff dimension dimj^(£J) and the random Hausdorff dimension dim|f (E) are linked by 
the relation 

KPZ formula. Almost surely, we have 

(3.3) ( (^l) =iimH{E) 
where £ is the so-called structure exponent of M: 

£(g) = (2 + ^)<?-^ 2 . 

Actually, Formula (3.3) remains valid for any MRM as constructed in (2.1) as soon 
as M possesses moments of negative order (see [7]). 

A further understanding of the geometry of the measure M is to define a Riemannian 
metric such that the corresponding volume form coincides with M. Such a structure 
permits to define fundamental objects such as distance, arclength, geodesies that are 
naturally associated to the measure M. This is an open problem suggested by K. Izyurov 
at the conference "Les Diablerets". 

Theorem 2.6 allows to understand the measure M on Br as the volume form of a 
Riemannian manifold (up to a set of null M-measure) . The geodesies are easily described 
via the transport maps. By sticking to the notations of Theorem 2.5, we define T = 
and x = T -1 - Since T is an isometry of metric spaces, a curve 7 : [0, 1] — > C is a geodesic 
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on C if and only if IX7) C Br is a geodesic for the Euclidian metric, that is a segment. 
The geodesic joining x, y € B C Br (parameterized by constant speed) is thus given by 

rf* : t G [0, 1] X (tr(x) + (1 - i)r(y)) € C. 

In what follows, we give a formal interpretation of the KPZ formula as well as a 
conjecture on the Hausdorff dimension of the geodesies. Consider a compact set K, 
contained in Supp(x)- Then x(K) is a subset of C but may also be thought of as a 
subset of Br. Thus we can define the Hausdorff dimension of x(K) as a subset of C 
equipped with the distance d, denoted by dim.j I (x(K)) and the Hausdorff dimension of 
x{K) as a subset of Br equipped with the Euclidian distance, denoted by dim|^(x(i^)). 
By isometry between (C,d,M) and (Br, d e , M(Br)Xr), we have the obvious relation 

dim d H { X (K)) = dim H (K). 

Furthermore, if we assume that formula (3.3) still holds for these Hausdorff dimensions 
(defined in terms of distances not measures), we will have 

The KPZ formula can thus be reinterpreted in terms of the graph of the function x : X 
sends a set K with Hausdorff dimension q towards the set x{K) with Hausdorff dimension 
£(g/2). In particular, the geodesies with respect to the metric d should have Hausdorff 

2 

dimension £(1/2) = 1 + 

We cannot give a rigorous argument to that conjecture because the general theory 
of optimal transport suffers a definite lack of strong estimates concerning the optimal 
transport maps. So we cannot prove the KPZ formula connecting the Hausdorff dimen- 
sions defined in terms of distances (not measures) and we leave that point as an open 
question. 



3.2 Multifractal random change of time 

Another important motivation for constructing an optimal transportation maps associ- 
ated with MRMs is the possibility of defining multifractal random change of times as 
suggested by Mandelbrot in the Id case: if B is a 1-d Brownian motion and M is a 1-d 
MRM, we can define the process t t— > #M([o,t]) • This process has been widely used in 
financial modelling since it enjoys nice properties: it is a square-integrable martingale 
with long-range correlated stationary increments and power law spectrum. Extending 
that construction may be done via optimal transportation maps. 

Let us fix R > and a m-dimensional MRM as constructed in 2.1. Let T,x be the 
maps pushing M forward to A/ ^P fl ^ dx and vice versa (sticking to the previous notations, 
we denote by B the support of T). We further consider a Gaussian white noise W(dx) 
defined on R m . Given x € M. m , we denote by C(x) the cube [0, x{\ x • • • x [0, x m ] (with 
the convention that, if Xi < 0, the interval [0,Xj] stands for [xj,0]). Finally, we define the 
process 

Vx G B R , B{x) = W(T(B n C{x))) . 
It is readily seen that the process B is isotrop and stochastically scale invariant: 

VA<1, (B(Xx)) xeBR ^X^e^(B(x)) x£B , 
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(a) Conformal geodesic ■y 2 ~ 2 




(b) Comparison of conformal/optimal transport geodesies for -y 2 ~ 1.5 

Figure 2: Simulations of conformal geodesies (weighted by the density of the MRM along the 
path) and geodesies obtained via optimal transport for a log-normal MRM. 
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where Q\ is an infinitely divisible random variable, independent of (S(x)) a;gS , the law 
of which is characterized by: 

E[e iqQx ] = \-^ q l 

Of course, we can generalize this approach to make a multifractal random changes 
of time to many other multidimensional stochastic processes like Levy noise, fractional 
Brownian motion, and so on... 

4 Proofs 

We first prove the key estimate of the paper. 
4.1 Fundamental result 

Let B be a Borelian subset of Br and k be a probability measure on Br supported by 
B, meaning k(B) = 1. We assume that the set B is equipped with a distance d and that 
the completion of B with respect to the distance d, denoted by (C,d), is compact. We 
assume that the Borelian subsets of B with respect to the Euclidian topology coincide 
with the Borelian subsets of B w.r.t. the distance d so that we can extend the measure 
k to the whole Borelian sigma-algebra of C by prescribing 

VicC Borelian, k(A) = k(A n B). 

The classes R Q ,R~. For any a > 0, we introduce the set R a of Radon measures v 
on C satisfying: for any e > 0, there are 5 > 0, D > and a d-compact subset K t C C 
such that v(C \ K e ) < e and the measure v t (dx) = lK e (x)v(dx) satisfies 

(4.1) Wd-open ball, v € {U) < D x Di&m d (U) a+s , 

where Diam^ denotes the diameter of the ball U with respect to the distance d. We 
further define the set of Radon measures R~ = C\p <a Rp- 

We further give an energy condition for a Radon measure v to be in the class R~. 
For a > 0, we define 

(4.2) CaM= JJ c _±_, {dx y {dy) . 
It is plain to see that 

v satisfies C a {y) < +oo => v € R~. 

Conversely a measure v obeying (4.1) satisfies Cp{v) < +oo for each /3 < a + 5. 

Notations: In what follows, we use the superscript e (i.e. R e ~ and CI) to mention 
that the distance d is equal to the Euclidian distance. □ 

Proposition 4.1. Assume that the measure k belongs to R e a . Let N be the Radon 
measure on Br defined, F-almost surely, as the following limit (in the sense of weak limit 
of Radon measures): 

N(dx) = ]xmN l (dx), where N^dx) d = e Ul{x) n(dx) . 

Z->0 

If ip(2) < a then the martingale (N^B)^ is uniformly integrable and, consequently, N 
is non-trivial and almost surely supported by B. 
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Proof. We remind the reader of the fact that the family (N^B)^ is a right-continuous 
positive martingale and thus converges almost surely. 

Since n belongs to R a , for each e > 0, we can find 5 > 0, D > and a compact subset 
K such that k(C \ K) < e and 

nmsup < —a — o uniformly for x G K, 

r—^oo V 

where B* stands for the d-ball of radius e~ r and center x. This implies Cp{Kx) < +oo 
for each a < /5 < a + S. 
Then we compute 



E 

where K\ is given on R m by 



Ni{BHK) 2 = / e^ K ^ y - x ^ K (dx)K(dy) 

- 1 JBnK Jbhk 



Ki(x) = [ pi(gx)H(dg) 
JG 



where pi(y) (y G R) is defined by ln(T/|y|) if I < \y\ < T, Pl (y) = ln(T/Z) + 1 - \y\/l if 



< I and otherwise. In particular, on B, K[ is not greater than C + In for some 
positive constant C. As a consequence, for some positive constant C which may change 
along the inequalities, we have 



supE 



NtiBnK) 2 } < C I [ e min ^K(dx)K(dy) < C'C^k) < +oo 
- 1 Jbdk Jbhk 



since i/j(2) < f3. The martingale (vi(B n K)^ is bounded in L 2 (Q) and is therefore 
uniformly integrable. It is plain to deduce that the martingale (z/j (£?)), is uniformly 
integrable. □ 

Theorem 4.2. 1) Assume that the measure k belongs to R a n R^ for some ^ > Y>(2). 
Then 

N g R aS ^m +fi2) _^ {iy 

Consequently, we also have: 

k G R~n R e ~ for some q > ip(2) =>- N G fl~ c _^ r21 

i?j In particular, we have in the Euclidian case: if k € i?^ (resp. R^~) and ip{2) < a 
tfien M G (neap. R^,^)- 

Proof. Since k belongs to i? a n R^, for each e > 0, we can find 5 > 0, D > and a 
compact subset -fT' for the d-topology (<i-compact for short) such that k(C \ K') < e/2 
and 

(4.3) limsup — — < —a — 5 uniformly for x G K', 

r— ¥oo T 

where B* still stands for the <i-ball of radius e~ r . 



13 



Since k belongs to i?|, for each e > 0, we can find a compact subset K" for the 
Euclidian topology such that k(Br \ K") < e/4 and Cp{l.K"{x)n{dx)) < +00 for each 
p < q. K" n B is a Borelian subset of B (for the Euclidian topology) and hence a 
Borelian subset of C. Since M is inner regular on C (recall that C is compact), we can 
find a (i-compact set K of C such that K C if" and M(K" \ K ) < e/4. Finally we set 
K = K ' n K, which is d-compact. The set K also satisfies: M(C \ K) < e, (4.3) is valid 
on K and (71 (Ismc) < +°o for any /3 < 

Even if it means multiplying k by a constant, we assume k(K) = 1. We consider on 
x K the probability measure Q defined by 



/(w,x)dQ = E / f{uj,x)N(dx) 

for all measurable nonnegative functions /. 
Given V < I, we define 

UJl'^x) = uji(x) - oj v (x). 

For a sequence £1 < • • • < Z n , the random variables . . . ,0Ji n _ 1) i n are Q-independent 

and 

e Xui i' l i dQ = K[e {1+X)ull » l i ix) ] = e ^(i+A)lnfeA)_ 

The process u € M+ i-> w e - u is therefore an integrable Levy process (we can consider a 
version that is right-continuous with left limits). From the strong law of large numbers, 
we have 



This implies that, P a.s. 



a.s., > tp (1) as u — > 00. 

u 



Jy a.s., > tp (1) as u — > 00. 



Therefore, P a.s., for each e > we can find a compact K\ C K such that N(K\K^) < e 
and ^e-^) _^ ^'(1) uniformly w.r.t. x € if* as u — >• 00. Now we define 

iVg(cfy) = lime a V- 9(2/) rc(dy) and PJx) = I NJdy). 

J B%C\K 

We further define the function g by 

\ otherwise. 
Thus we have P q {x) = f K 6 q (x,y)N q (dy) and 

P q dQ = E / / 0,(x, y)N q (dy)N(dx) 
Jk Jk 



< 



BnK J BnK 



(x,y)e m ( C+H tt ) ) K{dx)K{dy) 



<e^ c [ [ 9 q (x,y)e-^- 1 K (dx) K (dy). 

J BnK J BnK \y ~ x\ n} 



14 



By using the above relation, we obtain 

[j2^ n PndQ = J2 [ I e^-^ n Xn (y,x)- 1 K (dx)K(dy). 

J ^.WW \y-x\W) 

Note that (for some positive constant D) 

n>l l<n<-lnd(x,y) 1 ' 

in such a way that we obtain 

/ ^ n Pn dQ<D f [ 1 1 K(dx)K(dy) 

J ^ JBnKJBnK \y-x\n 2 > d(x,y)P 

We want to prove that the latter integral is finite for a well chosen (3. We fix 

if>(2) <(3 = (a + 5/2 f~^ 2) + ^(2). 
We consider p,q > 1 satisfying the relation - + - = 1 and given by 

<° = W) and q = 7^WY 

By Holder's inequality, we have 



< 



1 



1/p 



K(dx)n(dy) I I — — 7-5 — -t-^t- K(dx)n(dy) 

BnKJBnK \y-x\P*m v ;i y '\ lJ BnK J BnK d(x,y)i(^( 2 )) V 1 y U> 



1 



1/9 



< (^C ? e (l^nB(x)K( ( ix))J [C a+< 5 /2 (%nB(a;)K(^)) ; 

in such a way that the above integrals are finite. 

We deduce that, Q a.s., e^ n P n -> as n — > oo. Therefore, P a.s., eP n P n — s> as 
n — >■ oo iV-almost surely. So we can find a compact C C such that N(C \ K^) < e 
and 

111 t^fi \ X J „ _ if _ _ o 

hm sup < — p uniformly for x € iv e . 

Finally we can set K = Kj n K 2 and N^(dx) = Ig-(x)iV(<ix). We obtain 

lnjV*(fl£) In/x- n B ;^-»(%„((fe) 
hm sup = hm sup - 

n— >oo Tl n—too Tl 

<-/3 + V'(l) 

uniformly w.r.t. t & K. We have proved N £ | P-almost surely. □ 
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4.2 Composition of MRM 

Now we prove that a non-degenerate log-infinitely divisible random measure M, i.e. 
satisfying Y>(2) < +00 and < d, can be decomposed as an iterated multiplicative 

chaos. 

Since ip(2) < +00 we can find an integer n such that 
(4.4) mip(2) < n(m-ip'{l)). 

Then we can find n independent independently scattered log-infinitely divisible random 
measures . . . , associated to ((p/n,9) (remind of the definition in subsection 
2.1). We assume that the random measures ■ ■ ■ , are constructed on the prob- 
ability spaces (flWjPW ),..., (ft( n \p( n )). We define ft = f^ 1 ) x ••• x f]W equipped 
with the product <r-algebra and the product probability measure P = P^ 1 ) <g> • • • <g> . 
The corresponding processes u/j associated to fi^ , . . . , are respectively denoted by 



(i) , » 



, ujf 1 ' . Finally, we denote by E^ the conditional expectation given the variables 

We define recursively for k < n: 

M^{dx) = dx and M^(dx) = lim M^ k ~ l \dx), 

where the limits have to be understood in the sense of weak convergence of Radon 
measures. Note that the choice of n makes valid the relation 

mib(2) k , , . 

VJfc<n-l, < m- 

n n 

Hence, we can apply recursively Proposition 4.1 and Theorem 4.2 (with the distance d 
equal to the Euclidian distance and k G is the Lebesgue measure) to prove that, for 
each k < n, 

M(k) e ^Ivd) and Eik) l Mik) (B R )} = M^HBr) P a.s.. 

Thus we have AfW £ R e ~_^ {1) and E[M^(B R )] = X(B R ) (the Lebesgue measure of 
Br). 

What we now want to prove is that the measure has the same law as the measure 

Af(dx) = lim e u i 1){x)+ - +u i n){x) dx. 
i— >o 

We consider on Q the cr-algebra Qi generated by {u4 (x) , . . . , oj^ {x);x S P m , T > r > I}. 
The conditional expectation of M^ n '(E) w.r.t. Gi is easily computed since, for each k < n, 
the martingale (Afj (A))i is P( fc ) -uniformly integrable. Indeed, we have: 

E[MW(4)|g,] = E[E[MW(4)|^,. . ,^1 (^(i))^,^] \Qi) 
= E [ E (n) [ M ^(A)\(4 n \x)) xeMm>T>r>l ] \g{\ 

,(») 



E[ / e"i {x) M^-^{dx)\gi] 



A 



A 
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This latter quantity has the same law as Mi(A). Since the martingale (E[Af ( - n ^(A)|ft]) / is 
uniformly integrable, we deduce that the family {Mi(A))i is uniformly integrable. Hence, 
both random variables M(A) and M^ n \A) have the same law. In particular, M £ 

Corollary 4.3. If ij)(2) < +00 and < m, then M belongs to R 



e— 

m— 0'(1)' 



Remark 4.4. The same composition argument shows that if a measure k € R e a and if 
M is defined as the limit 

M(dx) = lim e Wl{x) K{dx) 

and satisfies ip(2) < +00 then M £ Though we won't use it in this paper, 

this is an important result concerning the structure of the support of multifractal random 
measures. 

4.3 Proof of Theorem 2.3 

In the case where ip{2) < +00 and ^'(l) < 1, we will prove that the measure M does 
not give mass to small sets. Let S C Br be a small set, that is a set with Hausdorff 
dimension (w.r.t. the Euclidian distance) not larger than m — 1. 

From Corollary 4.3, M belongs to the class fi^l^m- Since < 1, M then 

belongs to the class for some j3 > 0. We fix e > 0. So, P a.s., we can find a 

compact set K C Br and 5, D > such that M{Br \ K) < e and for all open balls 
UCB R : 

M(U n K) < Cdiam e ([/) m - 1+/3+<5 . 
Since m — 1 + /3 + S > dim// (5), we can find a covering of S by open balls (Ui)i such that 



m-l+/3+S < £ 



y^diam e ([/i)'' 
Then we have 

M(S) < M{S \K) + M{S n K) < e + ^ M(Ui n K) < e + ^ diam e ([/i) 1+/3+7 < 2e. 
As we can make e as small as we please, we deduce M(S) = 0. □ 
4.4 Proof of Theorems 2.5 and 2.6 

We proceed recursively to prove the existence of an optimal transport between the mea- 
sure Aft* -1 ) and Af( fc ) on some appropriate Riemannian manifold: 

1) Step 1: We focus on the measure M^ l \ which has structure exponent C (?) = 
mq — Relation (4.4) and the convexity of ip imply the following inequalities 

m < m _ L rnm _ 1 
n n m n m 



Hence we can apply Theorem 2.3 to find two optimal transport maps x^; 1 ^ 1 ) that 

to M (1) and vice versa. Furthei 

inf [ \T(x) - x\ 2 M (1) (dx) 



respectively push Xr forward to and vice versa. Furthermore, the quantity 



T:Br^b r J b 

T # M (1) =\ R 
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is achieved at r^. We can also apply Theorem 2.4 to find a compact Riemannian 
manifold (C^\ g^) and a Borelian subset of Br such that: 

-B^ is dense in Br for the Euclidian distance and has full M-measure, that is 

mW(b r \bW) = o, 

-C^ 1 ) is the completion of B^ with respect to the geodesic distance on 

-the volume form on coincides with the measure on £?' ) 

-in a system of local coordinates, the Riemannian metric tensor on reads 

gV=e^(u)(dxl + --- + dx 2 rn ) with e^(co) = M(l) ^ R)2 . 

Furthermore, from Proposition 4.1 and Theorem 4.2, € ^2~ip'(i)/n' This ends 

up the first step of the induction. 

2) Step 2: we assume that, for some k < n, we may find a compact Riemannian 

manifold {C^ k \g^ k ') and a Borelian subset B^ of Br such that: 

-B^ is dense in Br for the Euclidian distance and has full M-measure, that is 

M^{B R \B^) = 0, 

-C( fc ) is the completion of B^ k > with respect to the geodesic distance on C^ k \ 

-the volume form on coincides with the measure M( fc ) on B^ 

-in a system of local coordinates, the Riemannian metric tensor on reads 

g w = e^\u)(dxl + ■■■ + dx 2 j with e^ k \io) = M( ^ ( f fl)2 . 

We denote by = rW o ■ ■ -oTW : (C^ k \g^) — > (Br, d e ) the isometry constructed 
recursively. From Proposition 4.1, the measure M^ k+1 ^ is supported by B^ almost 
surely, that is M^ +1 \Br\B^) = almost surely, so that M( k+1 ^ extends to a measure 
on by prescribing: 

V.4 C C (fc) Borelian, M^ k+1 \A) = M {k+l) {A n B {k) ). 

Furthermore, from subsection 4.2, we have M<® 6 R e ~ k We now apply Theorem 

m-^V'(l) 

4.2 where the distance d is equal to the geodesic distance on C^ k \ denoted by d^ (the 
corresponding class R a will be denoted by Ra )■ Since is the volume form on 

we have 

M( k \u)=Dr m , for any open ball U with radius r. 

Hence G r!$ . The assumptions of Theorem 4.2 are thus satisfies with a = 

m — ^ip'(l) an d <; = m — ^V'(l) (and thus we have ip(2)/n < <j because of (4.4)). It 
follows that M( fc+1 ) € R^~ ^ {2) v>(a)-v'(i) ' because of (4.4) again, we have 

m ( 1_ mn-fc^'(l) ) H n 

V>(2) ^(2)-^(l) ^(2) , 

m(l ) H > m(l 



van — ki/}'(l) n mn — kip'(l) 

mip(2) 



> m 



> m — 1, 
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so that we can show, as in the proof of Theorem 2.3, that M^ k+1 ^ does not charge 
the small sets of {C^ k \g^). Hence we can apply Theorem A. 3 to find two optimal 

transport maps a( k+1 \ f3^ k+1 ^ that respectively push forward to M^ k+1 ^ and vice 

versa. Furthermore, f3^ k+l ^ can be rewritten as (ipW)- 1 o r( fc+1 ) o ip( k ) where : 
Br — > Br is the gradient of some convex function. The function p^ k+l "> = tp^ o /3( fc+1 ) = 
p( fc +!) o • • • o I^ 1 ) thus pushes M^ k+l ^ forward to Xr. Besides the quantity 



inf / S k \T{x),x) 2 M [k+l \dx) 

T:B R ^B R J B 



is achieved at f3^ k+l \ Since S k \x,y) = d e (p>^ (x), p^ (y)), we deduce that 
realizes the above infimum if and only if 

^(fc+i) = ^k) 

o /3( fc+1 ) realizes the infimum 
inf / \T{x)-^ k \x)\ 2 M {k+1 \dx). 



T # M {k+1) =X R 

Then we can apply the same machinery as in subsection 2.2 to construct the k + 1-th 
Riemannian structure, which ends up the induction. □ 



A Background about optimal transport theory 
A.l Monge problem 

We remind the reader of the following classical results 

Definition A.l. Push-forward of measures. Let [i,v be two measures respectively 
defined on the measured spaces E and F. We will say that a measurable mapping p> : 
F — > E pushes the measure v forward to n if both measures and v o p" 1 coincide. In 
that case, we write p#v = fi. 

Definition A. 2. Small sets. Given a metric space (X, d) with Hausdorff dimension n, 
a small set is a set with Hausdorff dimension not greater than n — 1 . 

Given two probability measures [i and v on Br, a coupling of (fi, v) is a probability 
measure tt on Br x Br with marginals [i and v. A coupling tt is said to be deterministic 
if there is a measurable map T : Br — > Br such that the map x G Br h-» (x, T{x)) pushes 
fx forward to tt. In particular, for all i/-integrable function ip, one has 

/ l p(y)dv(y)= ip(T(x))dn(x). 
Jb r Jb r 

Such a map T is called a transport map between // and v. 

The Monge-Kantorovich problem on the ball Br can be formulated as follows. Given 
a cost function c defined on Br x Br, one looks for a coupling tt of (/i, v) that realizes 
the infimum 

C{n,v) = inf / c(x,y)dn(x,y) 

JB R xB R 
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where the infimum runs over all the coupling ir of (//, u). Such a coupling is called optimal 
transference plan. If the coupling ir is deterministic, the corresponding transport map T 
is called optimal transport map. The optimal transport cost is then the value 

/ c(x,T(x))dfj,(x). 
Jb r xB r 

The search of deterministic optimal transference plans is called the Monge problem. 
A. 2 Solution to the Monge problem 

We have (see [10, Theorem 10.28] and [11, Theorem 2.12 iv] for the last statement) 

Theorem A. 3. Let X be a Riemannian manifold isometric (as a smooth Riemannian 
manifold) to the closed ball Br. We denote by f : X — > Br the corresponding isometry 
of Riemannian structures. Let c : X x X ^-M. be the cost function given by 

c(x,y) = d(x,y) 2 

and fi,v two probability measures on X. Assume that the measure /i does not give mass 
to small sets. Then: 

1) there is a unique (in law) optimal coupling tt of (fJ,,u) and it is deterministic. 

2) there is a unique optimal transport map T (i.e. uniquely determined \x almost ev- 
erywhere) solving the Monge problem. Furthermore, we can find a lower semi- continuous 
convex function <fi defined on Br such that 

T(x) = f 1 oV^o/(x) 

for every x G f~ 1 ({y £ K; </> is differentiable at y}). 

3) Supp(i/) = T(Supp(/x)). 

4 ) Finally, if v does not give mass to small sets either, then there is also a unique 
optimal transport map T 1 solving the Monge problem (of pushing v forward to fi). We 
can also find a lower semi- continuous convex function tp defined on Br such that 

T'(x) = f- 1 oVcf>of(x) 

for every x € f^ 1 [{y £M;^ is differentiable at y}) . T and T' satisfy, for pi almost every 
x € X and v almost every y € X , 

T'oT(x)=x, ToT'(y)=y. 

Proof. There is an easy way to deduce the above theorem from [10, Theorem 10.28]. 
Because of the isometry with the closed ball Br, the above theorem is basically of Eu- 
clidian nature. Indeed, it is plain to see that tt is an optimal coupling of fx, v for the cost 
function c(x,y) = d(x,y) 2 on X x X if and only if II = 7T#(/, /) is a coupling of the 
probability measures of ff^v on Br. In the same way, T : X X is an optimal 

transport map such that T#fi = v if and only if 9 = f o T o f" 1 is an optimal transport 
map such that 0^(f^^) = f^v for the Euclidian quadratic cost on Br. So the proof 
of the above theorem boils down to applying [10, Theorem 10.28] in the Euclidian case 
with quadratic cost function. It is then plain to complete the proof. □ 

Remark A. 4. In case 4) is satisified, it is more convenient to restrict the support ofT 
and T' respectively to {x £ X ; <fi is differentiable at f(x) and T' o T{x) = x} and {x G 
X;ip is differentiable at f(x) and ToT'(x) = x}. In that way, T : Supp(T) — > Supp(T') 
and T' : Supp(T') — » Supp(T) are both bisections. 
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